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The  Topside  Ionosphere  Plasma 
Monitor  (SSIE)  for  the 
Block  5D/ Flight  2 DMSP  Satellite 


1.  INTRODUCTION 

This  report  describes  the  Topside  Ionosphere  Plasma  Monitor  flown  on  the 
Block  5D/Flight  2 Satellite  of  the  USAF  Defense  Meteorological  Satellite  Program 
(DMSP).  The  plasma  monitor,  designated  as  Special  Sensor  Ion/Electron  (SSIE), 
consists  of  an  electron  sensor  and  an  ion  sensor  mounted  on  a boom  from  the 
DMSP  satellite  as  shown  in  Figure  1.  The  electron  sensor  measures  the  ambient 
electron  density  and  temperature  and  the  electrostatic  potential  of  the  vehicle  with 
respect  to  the  ambient  plasma.  The  ion  sensor  measures  the  densities  of  the 
major  ion  species  present,  the  average  ion  temperature,  the  average  ion  mass, 
and  the  vehicle  potential.  The  plasma  scale  height  at  the  satellite  is  determined 
with  the  use  of  both  instruments.  From  the  continuous  readout  of  data  from  the 
plasma  monitor,  the  plasma  scale  height  and  density  are  determined  as  a function 
of  time,  and  plasma  irregularity  spectra  can  be  calculated.  The  vehicle  potentials 
and  plasma  densities  obtained  independently  from  the  ion  and  electron  sensors  are 
compared  as  a check  of  instrument  performance  and  of  data  reduction  procedures. 

The  electron  sensor  of  the  SSIE  system  is  a spherical  Langmuir  probe.  The 
electron  sensor,  together  with  its  electronics,  is  shown  schematically  in  Figure  2. 
The  flight  hardware  and  the  electron  data  analysis  procedure  are  described  in 
Section  2.  The  ion  sensor  of  the  SSIE  system  is  a planar  electrostatic  particle 

(Received  for  publication  17  March  1978) 
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trap.  The  ion  sensor  and  its  electronics  are  shown  schematically  in  Figure  3. 

Tlie  flight  hardware  and  the  ion  data  analysis  procedure  are  described  in  Section  3, 
the  method  of  determining  the  plasma  scale  height,  in  Section  4,  and  the  SSIE 
telemetry  allocations  and  electronics,  in  Section  5.  The  time  sequence  of  events, 
bias,  sweep,  and  internal  calibrations  for  both  sensors  are  given  in  Table  1. 


ELECTRON  SENSOR 


AR  CELLS 


VERTICALLY  TO 
TIMES 

» PARALLEL  WITH 


ENSOR  PARALLEL 


LLEL  WITH  -X  AXIS 


Figure  1.  F2  Spacecraft  Showing  Location  of  SSIE  Sensors  (not  to  scale) 


t 


2.  THE  ELECTRON  SENSOR 


The  electron  sensor  consists  of  a 1.  75  in.  diam  spherical  gold  plated  aluminum 
collector  surrounded  by  a 2.  25  in.  diam  spherical  grid  made  of  gold  plated  tungsten 
wire  of  0.  800  transparency.  It  is  mounted  at  the  end  of  a 2.  5— ft  rigid  boom  deployed 
after  the  spacecraft  is  stabilized  in  orbit.  A cross  section  of  the  electron  sensor 
with  its  mounting  assembly  is  shown  in  Figure  4.  The  electron  and  ion  sensors 
attached  to  the  boom  are  shown  in  Figure  5. 

The  electron  sensor  has  two  modes  of  operation:  In  Mode  1,  a fixed  voltage 
with  respect  to  spacecraft  ground  is  maintained  on  the  outer  grid  in  order  to  obtain 
continuous  electron  density  measurements.  The  fixed  voltage  is  either  +1.  5 V 
(Bias  Mode  1 in  Figure  2)  or  0 V (Bias  Mode  2)  with  respect  to  spacecraft  ground, 
and  it  may  be  changed  by  ground  command.  In  Mode  2,  a linear  sweep  voltage 
with  respect  to  spacecraft  ground  is  applied  to  the  outer  grid  as  shown  in  Figure  2. 
This  mode  gives  electron  temperature  and  the  vehicle  potential  by  a Langmuir 
probe  type  analysis  of  the  current  collected  versus  applied  voltage.  At  all  times, 
the  voltage  on  the  collector  is  +20  V with  respect  to  the  outer  grid  to  ensure  that 
all  electrons  passing  through  the  outer  grid  will  be  collected.  Mode  1 and  Mode  2 
operations  of  the  electron  sensor  are  automatically  alternated  by  the  SSIE  Time 
Sequencer  as  described  in  Table  1 and  Section  5.  1.  One  Mode  2 sweep  operation 
of  10-sec  duration  is  performed  every  128  sec,  and  Mode  1 operation  is  maintained 
during  the  remainder  of  the  128-sec  cycle.  The  values  for  the  applied  voltages 
during  the  electron  sweep  mode  are  given  in  Table  2. 

2.1  Electron  Sensor  Data  Analysis 

Analysis  of  the  electron  sensor  data  starts  with  the  determination  of  the  elec- 
tron temperature  Te,  the  vehicle  potential  <j>g,  and  the  electron  number  density  Ng 
from  analysis  of  the  Mode  2 sweep  data.  It  is  assumed  that  Tg  and  <j>  are  constant 
for  the  next  128  sec  and  they  are  used  to  determine  the  instantaneous  values  of  the 
ambient  electron  density  from  the  current  flowing  to  the  electron  sensor  in  the 
fixed  voltage  mode  (Mode  1). 

2.2  Sweep  Mode  Analysis  (Mode  2) 

The  current  (Ie>  to  a spherical  electron  sensor  is  a function  of  the  sensor 
potential  (b)  with  respect  to  the  plasma.  If  <t>  is  positive,  then  electrons  are  accel- 
erated to  the  sensor;  if  <t>  is  negative,  they  are  retarded  as  they  approach  the  sen- 
sor. The  expressions  for  the  current  reaching  the  sensor  and  its  first  derivative 
are  given  here,  assuming  a Maxwellian  distribution  function  for  the  ambient  elec- 
tron population  for  a retarding  potential  (b  < 0): 
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Table  1.  Time  Sequence  of  Events,  Bias  and  Sweep  Voltages,  DMSP,  Elock  5D,  F2 




Event 

Time  (sec) 

Event 

Applied  Signal 

Monitor 

Level 

Duration 

(sec) 

First 

Repeats 

Every 

Calibration  1 

+5  X 10'8a  (Ion) 

-1  X 10_6a  (Elec) 

0.  2 V 

2 

0-2 

1024 

Calibration  2 

+5  X 10'Ua  (Ion) 

-1  X 10'9a  (Elec) 

0.  4 V 

2 

2-4 

1024 

Electron  Sweep 
Flag 

None 

0.  6 V 

2 

8-10 

128 

Electron  Sweep 
(Outer  Grid) 

+5  V -»  -8  V 

+5  V - +1  V 

10 

10-20 

128 

Ion  Sweep 

Flag 

None 

0.  8 V 

2 

24-26 

128 

Ion  Sweep 
(Outer  Grid) 

-5  V — +12  V 

+1  V - +5  V 

12 

26-38 

128 

6 - potential  of  sensor  with  respect  to  plasma  - (volts) 

s p 

(Js  = potential  of  vehicle  with  respect  to  plasma  (volts) 

<t>p  = potential  of  sensor  with  respect  to  vehicle  (volts) 

2 2 

A = surface  area  of  sensor  (m  ) = 4*r  , where  r is  sensor  radius  = 
0.  028575  m 

-19 

e = electron  charge  = 1.  602  X 10  coulomb 
a = sensor  transparency  = 0. 800 

_3 

Ne  = ambient  electron  density  (m  ) 

me  = electron  mass  = 9.  1085  X 10  ^ kg 

1/2  - 1 
a = (2  kTe/me)  ' = most  probable  electron  speed  (m  sec  ) 

Tg  = electron  temperature  (°K) 

k = Boltzmann  constant  = 1.38054  X 10  ^ joule/°K 


To  obtain  T . we  make  use  of  the  linear  relation  between  Log  (-1  ) and  d>  in 

0 0 P 

the  retarding  portion  of  the  I vs  <t>  curve.  (We  will  denote  LogjQ  by  "Log"  and  use 
’In1  for  natural  logs.)  By  substituting 


9Log(-ie)/a^p  = (ln(io))-1  a(-ie)/ie  d<t>p 


and  Eq.  (1)  into  Eq.  (2),  we  obtain 


3 Log  (-Ie)/34>  = e/k  Te  In  (10) 


T 5040  °K 
le  ' S 


where  S is  the  slope  of  the  semi-log  plot  of  I vs  ^ . 

6 p 

In  Figure  6,  Log  (-Ie)  and  its  derivative  are  plotted  as  a function  of  </>p.  The 
point  4>q  indicated  on  the  curve  is  the  point  at  which  the  probe  is  at  plasma  potential 
(<t>  = 0)  and  it  is  the  intersection  of  the  exponential  (retarding)  and  linear  (accelera- 
ting) portions  of  the  I vs  <t>  curve.  From  this  point,  we  obtain  the  value  of  vehicle 
e p 

potential  4>g  with  respect  to  the  plasma 


Figure  6.  Electron  Sensor  Theory:  Logjg(I)  vs  and  First  Derivative 


Using  the  current  1 measured  when  the  probe  is  at  plasma  potential  (<f>  = 0),  one 
calculates  the  value  of  the  ambient  electron  density  Ng  at  the  time  of  the  sweep 
from  Eq.  (1). 

N (t  ) = 2 I /Ae  a a . (10) 

e o o 

Equations  (1)  through  (4)  are  based  on  the  assumption  that  |ed|  <,kT.  In  fact, 
the  range  of  the  swept  voltages  may  exceed  the  range  of  voltage  where  this  assump- 
tion is  valid.  In  the  accelerating  potential  region,  the  finite  Debye  length  (length 
that  the  probe's  electrostatic  field  can  reach  into  the  plasma)  causes  the  current 
to  the  probe  to  be  limited.  The  slope  of  the  (-Ig)  vs  <£  curve  is  less  than  that  given 
by  Eq.  (4)  when  <j>  - $ <1.0  V.  In  the  retarding  potential  region,  the  measured 
current  to  the  probe  can  be  significantly  more  than  the  value  given  by  Eq.  (1)  when 
<t>  - <j>o  < -2  V,  due  to  high  energy  components  of  the  ambient  plasma,  secondary 
and  photoelectrons  from  the  grid  surfaces,  and  the  sensitivity  limits  of  the  instru- 
ment electronics.  The  data  analysis  procedure  in  Section  2.4  is  designed  to  avoid 
the  ranges  of  the  potential  for  which  the  foregoing  assumption  |e^|  < kT  is  invalid. 
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2.3  Fixed  Voltage  Mode  Analysis  (Mode  1) 

Once  and  Tg  have  been  determined  from  analysis  of  the  sweep  mode  data, 

one  can  proceed  to  derive  N at  all  other  times  if  it  is  assumed  that  T and  <4  do 

e e s 

not  change  during  the  128-sec  interval  between  sweeps.  During  this  interval,  the  * 

0 V or  +1.  5 V resting  bias  is  applied  to  the  grid  of  the  electron  sensor.  If  we 
designate  this  resting  bias  by  R volts,  then  the  electron  density  at  any  time  t dur- 
ing the  fixed  voltage  Mode  1 operation  is 

Ne(t)  = Ie(t)  Ne(to)/Ie(R)  (11) 

where  Ie(t)  is  the  sensor  current  at  time  t,  Ig(R)  is  the  sensor  current  during  the 
sweep  when  - R volts,  and  Ne(tQ)  the  electron  density  derived  from  the  sensor 
current  with  the  probe  at  plasma  potential  in  Eq.  (10). 

2.4  Electron  Data  Reduction  Procedure 

The  steps  followed  for  the  analysis  of  the  electron  data  are: 

1.  Take  the  Mode  2 measured  current  I as  a function  of  <f>  . From  the  start 
of  the  sweep  in  time  (d^  = +5  V),  fit  a straight  line  to  3 data  points  to  determine 
the  slope,  and  then  move  the  fitting  process  through  the  entire  sweep,  one  point  at 
a time.  (Measured  currents  corresponding  to  the  maximum  voltage  on  the  tele- 
metry channel  are  disregarded  because  the  true  current  is  out  of  range. ) We  find 
the  potential  <£LIN  at  which  the  absolute  value  for  the  slope  of  a 3 -point  fit  cen- 
tered on  iJ>LIN  is  a maximum.  If  two  or  more  fits  give  the  same  maximum,  choose 
the  fit  for  the  most  negative  value  of  <t>.  .... 

LjIIN 

2.  Take  Log  (-Ie>  as  a function  of  </>p  and  do  same  3-point  fit  as  in  step  1 to 
find  the  maximum  absolute  value  of  the  slope.  This  time  if  two  or  more  fits  give 
equal  maximum  slopes,  take  the  first  maximum.  Let  the  slope  of  this  line  be  S, 
at  a potential  of  <£. 

LaJVj 

3.  Compute  electron  temperature  Tg  from  Eq.  (8). 

4.  Compute 


^LIN  + ^LOG 

2 


5.  Interpolate  to  find: 

measured  current  I , at  6 - <4 

O P 0 

measured  current  IR,  at  = R 


(12) 
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6.  Compute  electron  density  Ng  from 


2 I /J  X 10”  „ 

N = -I- cm'3 

e Ae  a a 


(13) 


(T  obtained  from  step  3 is  used  to  compute  a. ) Numerically, 


N = -6.  896  X 109  I S cm'3 
e o 


(14) 


Thus  we  now  have  Ne(to>,  T0(t  ) and  tf>g(t  ) at  a time  t which  can  be  equated  to  the 
midpoint  of  the  sweep 

t + t. 


t = 
o 


s e 


(15) 


where  t and  t are  times  of  start  and  end  of  sweep,  respectively. 

7.  Calculate  the  electron  density  during  Mode  1 operation  using  Eq.  (11). 

The  density  is  calculated  during  the  following  118  sec  of  Mode  1 (M1EL) 
operation  by  assuming  that  the  temperature  Tg  and  vehicle  potential  <j>  are  constant 
between  sweeps  (Mode  2). 


3.  THE  ION  SENSOR 


The  ion  sensor  is  a planar  aperture,  planar  collector  sensor,  a cross  section 
of  which  is  shown  in  Figure  7.  It  is  mounted  on  the  2.  5-ft  rigid  boom  shown  in 
Figure  5 and  located  so  that  it  faces  into  the  spacecraft  velocity  vector  at  all  times. 
The  1.  00  in.  diam  aperture  is  covered  by  a double  grid  of  gold  plated  tungsten  wire. 
The  double  grid  is  used  to  minimize  electric  field  leakage  that  would  disturb  the 
trajectories  of  incoming  positive  ions.  Like  the  electron  sensor,  the  ion  sensor 
has  two  modes  of  operation  which  alternate  automatically  (see  Figure  3 and  Table  1). 
In  Mode  1,  the  retarding  (swept)  grids  in  Figure  7 are  electrically  connected  to  the 
spacecraft  frame  and  all  ions  passing  through  the  aperture  are  collected.  The  ion 
sensor  Mode  1 is  used  to  measure  instantaneous  ion  densities.  The  retarding  grids 
are  .--wept  once  every  128  sec  from  -5  V to  +12  V (Mode  2)  to  obtain  ion  tempera- 
ture. This  swept  grid  is  also  a double  grid  to  improve  the  uniformity  and  accuracy 
of  the  applied  potential  across  its  surface.  Between  the  retarding  grids  and  the 
collector,  a third  grid  is  maintained  negative  with  respect  to  the  collector  at  all 
times.  This  acts  as  a suppressor  to  the  photoelectrons  and  secondary  electrons 
generated  at  the  collector  by  solar  UV  and  energetic  particles,  respectively. 
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3.1  Sweep  Mode  Analysis  (Mode  2) 

The  current  flow  (I)  to  the  (planar)  ion  sensor  can  be  expressed  as  a function 
of  the  applied  potential  (4>p)  when  the  retarding  grid  potential  (<t>)  with  respect  to 
the  plasma  is  negative  (accelerating)  or  positive  (retarding),  assuming  that  the  out- 
ward normal  to  the  sensor  looks  along  the  spacecraft  velocity  vector.  This  analy- 
sis assumes  that  a single  Maxwellian  temperature  applies  to  both  ion  species. 

For  e<j>  >:  0 (retarding) 


~t  I N'  [ 


a.  exp  ( -xf ] 
1 + erf(x.)  + L 

1 Vs  ^ 


91 Ae2  a ^ Ni  2S 

■» T~  / exp  (-x.  ) 

% ,/ZWT  1 


where 


s » m. 

i 


erp  (-z  ) dz  (definition) 


| i>  = potential  of  sensor  with  respect  to  plasma  - <t>  + <t>  (volts) 

s p 

4S  - potential  of  vehicle  with  respect  to  plasma  (volts) 

<£p  = sensor  potential  with  respect  to  vehicle  (volts) 

2 2 

A = sensor  aperture  area  (m  ) = *r  , where  r is  aperture  radius  = 0.0127  m 

-1 Q 

e = electronic  charge  = 1.  602  X 10  coulomb 
a = sensor  transparency  = 0.  590 
N.  = ambient  ion  density  of  Ith  constituent  (m-3) 
m.  = mass  of  i1*1  constituent  (kg) 

a.  = (2kT/m.)1/2  = most  probable  speed  of  ith  constituent  (m  sec'1) 

T = average  ion  temperature  of  all  constituents  (°K) 
k = Boltzmann  constant  = 1.38054  X 10'23  joule/°K 
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j = number  of  ion  constituents 

Vg  = vehicles  speed  (m  sec'1)  = 7437  m sec'1 

For  e<t>  < 0 (accelerating)  on  the  retarding  grid,  the  current  collected  is  only  a 
function  of  the  flux  of  ions  through  the  aperture.  Put  4 = 0 in  retarding  Eq.  (16) 
and  x.  = V /a.  in  Eq.  (18).  Then  I is  constant  (not  a function  of  6 ) and  so 

1 S I Tp 


3.1.1  TWO  SPECIES  CASE 

To  handle  this  complicated  expression,  it  will  be  assumed  that  for  the  DMSP 
satellite  altitude  (circular  orbit)  of  835  km  we  have,  at  most,  two  ion  species, 
namely,  H and  0 . Of  course,  we  must  allow  for  either,  or  both  to  be  present; 
thus  for  j = 2,  nij  = 1.  66035  x 10'27  kg,  m2  = (16)  x (1.  66035  x 10'27  kg)  and  for 
j = 1,  either  ir^  = (1)  mp  or  = (16)  mp,  where  m = 1.  66035  x 10"27  kg. 

To  gain  an  understanding  of  the  behavior  of  the  current  flow  as  a function  of 
the  applied  voltage  on  the  retarding  grid,  let  us  for  the  moment  consider  what  the 
response  would  be  if  Vg  » a (very  low  ion  temperature  T),  and  we  have  two  com- 
ponents present,  where  m„  > m.. 

For  0 < <t>  < mjVg/2  e , and  xg  » 1,  erf  (Xj)  = erf  (x2>  = 1,  exp  (-x2)  = 
exp  (-x2)  = 0,  we  find  that 

I = Ae  Vg  o (N1  + N2) 

For  mi Vg/2  e < tj>  < m2V2/2  e,  Xj  « -1,  x2  » 1,  erf  (Xj) 
exp  (-x2)  = exp  (-x2)  = 0 we  find  that 

I = Ae  Vg  QN2  (20) 

For  <*  > m2V2/2  e , Xj  and  xg  « -1,  erf  (Xj)  = erf  (x2)  = -1,  exp  (-x2)  = 
exp  (-x2)'  = 0,  we  find  that 

I = 0 

These  three  operating  regimes  are  shown  graphically  in  Figure  8. 


(19) 

= -1,  erf  (x2)  = 1, 


17 


AeVsa  (Ni  + N2> 


Ae  Vs  <t  n2 


m*vs* 


Figure  8.  Ion  Sensor  Current  vs  Applied  Voltage  for  Two  Species 


Since  the  orbit  is  circular,  V is  constant  = 7.437  km  sec  (r,  orbital  period 

s 

= 101  min  28  sec),  and  we  can  calculate  from  Eq.  (18) 


= (1-  660  X 10~27)(7.  437  x 103)2  _ „ 287  y 
2x  (1.  602  x 10'19) 


<t>2  = (16)  ^ = 4.  586  V . 


These  are  the  two  points  that  give  x^  = 0 and  x2  = 0,  respectively.  It  can  easily 
be  seen  from  Eq.  (18)  that,  because  l/a[  occurs  as  a multiplier  on 
[Vg  -\/2  e^/m.],  and  02  give  zero  for  Xj  and  x2,  respectively,  for  all  at  (all  T). 
Examining  the  second  derivative  of  I (for  one  species),  we  note 


82I  _ -Ae2  <*N 


2 e x e~x 
*m  </>  * 


2 2 

it  is  obvious  that  8 I/8<A  = 0 when  x = 0;  hence  the  first  derivative  reaches  a 

P 

minimum  value  at  x = 0. 


! 
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Thus  as  the  temperature  is  increased  from  the  previously  assumed  small 
value,  the  slopes  at  ^ and  <\> 2 increase  from  -®,  but  the  slope  always  reaches 
its  maximum  negative  value  at  and  ^2  anc*  can  be  directly  related  to  the  density 

and  mass  of  each  species  and  their  temperatures.  For  illustrative  purposes,  Fig- 

2 2 3 -3 

ure  9 shows  a plot  of  I,  9I/90p  and  9 1/9 0pl  assuming  equal  densities  of  10  cm 

for  each  species  with  a temperature  of  3000°K  and  vehicle  potential  of  -2  V (the 

curves  were  computed  with  V = 7.  22  km  sec  1 as  opposed  to  correct  value  of 
-1  3 

7.437  km  sec  , but  this  merely  shifts  <K  and  </>„  a slight  amount  and  does  not 

* 1 “ 

change  the  character  of  the  curves).  It  can  be  seen  from  Figure  9 that  I is  deter- 
mined solely  by  species  (0+)  when  <j>  > 4.  5 V.  Because  of  the  behavior  of  erf  (x) 

2 P 

and  exp  (-x  ),  the  slope  9I/9ii  is  determined  solely  by  one  species  or  the  other  in 

P + + 

the  vicinity  of  <t>.  or  0„;(that  is,  by  H when  4 <4  V and  by  0 when  <4  > 4.  5 V). 
it  p p 

The  relative  contributions  of  the  two  species  can  be  clearly  separated. 

Figure  10  illustrates  Log  I vs  applied  voltage  for  the  same  data.  For  com- 
parison, the  electron  sensor  response  of  Log  (I)  as  a function  of  <4p  is  also  shown 
in  Figure  10.  It  is  clear  that  no  useful  data  can  be  extracted  for  the  ion  Log  I 
curve  (as  opposed  to  the  case  for  electrons  where  the  slope  is  inversely  propor- 
tional to  Te). 

Writing  out  the  separate  contributions  of  the  two  species  at  = 0 and  x2  = 0, 
respectively,  we  have: 

at  Xj  = 0,  <t>  = <4>j  + 4>g 


(24) 

(25) 

(26) 

(27) 


Curves  for 


the  single  species  case  are  denoted  by  (a). 
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Thus,  identifying  the  values  of  dp,  I and  91  9dp  at  which  the  slope  S is  a maximum 
negative  value,  we  obtain  dj,  Ij,  9Ij/9dp,  d.>,  I,  and  91,  9dp. 

From  Eqs.  (26)  and  (30),  we  obtain  the  vehicle  potential 


V 

s 


mp  (1  + 16)(dj 


>2>] 


(32) 


To  get  temperature  T and  the  densities,  we  shall  make  use  of  Eqs.  (28)  and  (29) 
only,  because,  as  can  be  seen  from  Figure  9,  the  first  derivative  at  dj  is  quite 
sharply  peaked  and  a small  error  in  location  dj  will  lead  to  a large  error  in  deduc- 
ing 9Ij/9dp.  At  d 2>  however,  the  peak  is  much  broader,  leading  to  a much  greater 
accuracy. 

Solving  Eqs.  (28)  and  (29)  for  N.,  and  T,  we  get 


N 


2 


16  m S„  it 

y- 


Ae  a 


1/2 

2I2e  1 

16  m S„  ir 
p 2 -J 


(33) 


and 


(Ae  aN2)2  e2 

2>rk  16  m S2 
P 


(34) 


For  Nj  we  could  use  Eq.  (24),  since  T is  now  known,  hence  aj  is  known,  but  again, 
since  the  slope  at  dj  is  so  steep,  the  value  of  Ij  will  be  susceptable  to  error.  A 
more  accurate  value  of  will  be  obtained  by  reading  off  the  current  Iq  at  dp  = -ds 
(we  now  know  the  value  of  ds),  that  is,  the  current  at  d “ 0 given  by 


Ae  V a 

3 


N, 


1 + erf  (V  /a. ) + 
s'  1 


2 2 
exp  (-Vg/a1) 

V 


+ Nr, 


1 + erf  (Vg/a2)  + 


a2  exp(-V2/a2) 
V 

s 


(35) 


where  aj  and  a2  are  defined  in  Eqs.  (27)  and  (31). 


4 


I 


1 


i 

i 


i 
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Thus 


2 I 

2-  . v 

Ae  V o 2 


+ erf  (V  a,,)  + 

S 1 


*9  exp  ( -V“  a~) 


v T 


1 + erf  (V  . a^)  + 


/ ,t2  9. 

a,  exp  ( - V ar) 

1 S 1 


(36) 


We  can  simplify  this  expression  slightly  by  noting  that  (Vg  a.,)  has  the  smallest 

value  for  largest  T.  If  we  take  T 103 4  °K  which  is  an  extremely  high  value,  then 

-1  “ “ 


a 2 = 3224  m see 
4.  884  X 10~'5  giving 


(V  a„) 

a ^ 


2.307,  erf(V./a9)  = 0.999,  and  exp  (-V“  a9) 


1 + erf  (Vg/a2)  + 


, ,r2  / 2. 

a9  exp  (-Vg/a2) 
Vs  ^ 


- 2.  000094  * 2 


For  a lower  temperature,  (Vg/a2)  becomes  larger,  giving  a value  even  closer  to  2. 
Because  (V  /a j)  = 0.  577  and  erf  (Vg/aj)  « 0.  585,  we  cannot  reduce  Eq.  (38)  any 
further,  and  thus 


2 I 


Ae  V 


2N, 

a 2 


1 + erf  (V  /a^ 


aj  exp  (-Vg/aj) 


(37) 


3.  1.  2 SINGLE  SPECIES  CASE 

The  case  where  just  one  ion  species  is  present  cannot  be  handled  by  the  pre- 
vious analysis  because  that  depends  on  locating  two  maximum  negative  slope  points 
where  the  first  .point  at  <t>.  can  be  identified  as  associated  with  H+,  and  the  second 
point  b2  with  0 . If,  however,  only  a single  minimum  slope  point  is  identified,  it 
is  not  known  a priori  whether  this  point  is  <j>^  or  b2.  To  illustrate,  in  Figure  10, 
is  we  identify  a maximum  negative  slope  only  at  </>  = 6.321  volts,  we  do  not  know 

+ P + 

whether  the  species  is  0 , giving  a vehicle  potential  dg  = -2  volts  or  H , giving  a 

vehicle  potential  of  -5.  052  volts.  The  corresponding  temperatures  work  out  to  be 

3000°K  and  1330°K,  respectively.  Essentially,  we  have  to  choose  a single  set  of 

equations,  (24),  (25),  (26)  or  (28),  (29),  (30),  depending  upon  whether  the  mass  is 

nip  or  16  nrip  (H+  or  0+).  The  clue  to  choosing  the  correct  set  is  to  make  use  of  the 

fact  that  when  <t>  < 0,  the  current  is  constant: 
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I 


A e V a 

I = § v 

o 2 


1 + erf  (V  a,)  + 

S ' 1 


a j exp  (-V^/aj) 


(38) 


Here  Nj  and  at  refer  to  the  now  single  species.  This  equation,  together  with 

equations  at  6 - <4,  (maximum  negative  slope)  = e> 

P 6 ^ wmm 


Ae  V oN 


V.  or 


91 
o 

d<j> 


-Ae  »N 


V 


t kT  m 


V 


W , + <t> ,) 

2 e — m 2.  o 


(39) 


(40) 


(41) 


are  solvable  numerically  for  the  four  unknowns  T,  nij,  and  <t>  , but  certainly 

not  algebraically  because  the  equation  in  I is  nonlinear. 

o 

There  are  two  ways  around  this  difficulty  and  both  ways  will  be  used  in  the  data 
reduction  scheme: 

1.  Make  use  of  Ng  derived  from  the  preceding  electron  sweep,  and  set  Nj  N 
in  Eq.  (38). 

2.  Assume  that  the  ion  temperature  T is  going  to  be  less  than  3000°K.  Then 
the  maximum  value  of  aj  will  be  given  by 


2k(300°!  . 7.063  * IQ3  , 

m * 


giving 


a exp  <-V?/a?) 
1 + erf  ( Va  /a. ) + — ? — i- 

S 1 ,,  r 


= 2.  0404  = 2 . 


This  is  a worst-case  and  value  gets  closer  to  2 for  T < 3000°K  and  is  very  close 

for  0+  (m,  = 16  m ). 

1 P 
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With  this  approximation,  Eq.  (38)  becomes 

I = Ae  V oN, 
o si 

and  is  independent  of  mass  giving  directly 

N,  = I /Ae  V cr 
1 o'  s 

Either  case  (1)  or  (2)  yields  the  value  of  N,  which  is  then  used  in  Eqs.  (39)  and 
(40)  to  give 


(Ae  aN^)  e 


ll  f Ae  V ttN,' 

2S1  - [h  — i—1 


kg 


(42) 


or,  in  atomic  mass  units 


M, 


(Ae  oN, ) e 


2SX  n 


h- 


Ae  V oN, 
s 1 


AMU 


(43) 


This  value  is  then  used  to  calculate 


T = 


(Ae  oNj)2  e2 
2~ 

2rrk  M.  m„  S. 

1 p 1 


(44) 


and 

(45) 

3.2  Ion  Data  Reduction  Procedure 

The  foregoing  analysis  then  leads  to  the  following  data  reduction  procedure: 

1.  Take  the  measured  current  I as  a function  of  <t>^.  Working  from  the  start 
of  the  sweep  in  time  (^  = -5  V)  and  disregarding  current  readings  when  the  tele- 
injtry  scale  is  saturated,  one  should  fit  a straight  line  to  3 data  points  to  determine 
the  slope  and  then  move  the  fitting  process  through  the  entire  sweep,  one  point  at  a 
time. 
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2T  M1  mp  - 


Find  the  center  point  at  which  the  slope  reaches  its  first  maximum  negative 
value  (<2^)  and  its  second  maximum  negative  value  (<2>2 ),  and  let  the  slopes  be  Sj 
and  S2>  respectively.  Take  also  the  currents  Ij  and  I2  at  <jj  and  d2-  If  only  one 
maximum  negative  slope  can  be  located,  proceed  to  step  4;  otherwise  go  to  step  2 

2.  Compute  vehicle  potential  <$  , given  by  Eq.  (32) 


17  V2  m 

2 e~P  ~ (*1  + *2>. 


and  read  off  (interpolate)  measured  current  (I  ) at  d -6 

ops 

3.  Compute  N2,  Nj,  T given  by  Eqs.  (33),  (37),  (34) 


N = 
2 


!?a 


S2  ” 1 

4 

m 1 
> 1 

V2  2 I e 

,3  + 2 

1/2 

* 1 

1 1 

4^  16  m S„  ir 

L p 2 

[•  2Iq  1 / [ a exp(-V2/a2) 

N!  =[a^-2N2J/  [l  + *r'tW+±-  -1  J 


(Ae  oN,)2  e2 

T = 2 

2 k k 1 6 m S, 
P & 


and  from  and  N2  we  obtain  the  average  ion  mass 


M_ 


(N.  + 16  N„) 

~^r 


AMU 


(46) 


Go  to  step  6. 

4.  Read  off  the  current  (I  ) at  <4  = 0 volts  and  compute  N,  = I /Ae  V . 

op  1 o'  s 

5.  Compute  Mj,  T and  <j>s  given  by  Eqs.  (43),  (44),  and  (45), 


M, 


2 Sj  ir  I 


(Ae  aNj)' 


Ae  Vs  rfO 


AMU 
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2(Ae  oN.  )2  e2 
- - i °j£ 

2 TrkM,  m S? 

1 p 1 


V?M1  »»„ 

5 1 -P  - <* 

2 e lj 


I i,  volts 


Repeat  step  5 with  Nj  = Ng  (value  from  previous  electron  sweep  M2EL)  generating 

a second  set  of  values  M^e,  Te,  The  two  sets  of  values  for  the  single  species 

ion  mass,  temperature,  and  vehicle  potential  can  be  compared.  Ideally,  the  two 

sets  will  be  identical  but,  in  fact,  they  occasionally  differ  and  both  sets  of  values 

are  made  available  for  a later  decision  as  to  the  best  values. 

6.  We  now  have  a set  of  ion  parameters.  Since  the  center  of  the  ion  Mode  2 

sweep  occurs  17  sec  after  the  electron  Mode  2 sweep,  these  ion  parameters  are 

assumed  to  be  valid  at  time  = t +17  sec.  Thus. 

o 

N.(t  + 17  sec)  = N,  + N„  or  N,  or  N 
i o l 2 1 e 

Mj(to  + 17  sec)  = Mgve  or  Mj  or  M ^ 

T.(t  + 17  sec)  = T or  Te 
1 o 

<t)  (t  +17  sec)  = <t>  or  <j>  e 
so  s S 

where  t is  the  time  of  the  center  of  the  previous  electron  sweep.  During  the  ion 
Mode  1 period,  the  instantaneous  total  ion  density  is  obtained  from  Eq.  (19)  as 

N(t)  =Nj(t)  + N2(t)  =I(t)/AeVga  (47) 


3.3  FJectrmtatic  Field  Penetration 

1 2 

It  has  been  shown  ' that  the  electrostatic  potential  in  the  spacings  of  a planar 
sensor's  grids  is  a function  both  of  the  electrostatic  potential  applied  to  the  grid 
wires  and  of  the  electrostatic  potential  of  adjacent  grids.  The  potential  in  the  spacing 


1.  Hanson,  W.  B. , Frame,  D.  R.,  and  Midgley,  J.  E.  (1972)  Errors  in  retarding 

potential  analyzers  caused  by  nonuniformity  of  the  grid-piane  potential, 

J.  Geophys.  Res.  77<No.  10):1914-1922. 

2.  Goldan,  P.  D. , Yadlowsky,  E.J.,  and  Whipple,  E.C.,  Jr.  (1973)  Errors  in  ion 

and  electron  temperature  measurements  due  to  grid  plane  potential  nonuni- 
formities in  retarding  potential  analyzers,  J.  Geophys.  Res.  78(No.  18): 
2907-2916.  
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H 

' 


is  the  effective  potential  of  the  grid;  it  should  be  used  in  the  above  data 

analysis  equations.  The  effective  potential  on  a grid  "i"  is  approximately  given  by 


Vi  + e (Vi-l 


-2Vi+Vi+l) 


where 


(48) 


V.  = the  applied  potential  to  grid  i 

the  applied  potential  to  the  grid  in  front  of  grid  i 
Vj+j  = the  applied  potential  to  the  grid  behind  grid  i 
e = field  penetration  factor 


In  theory,  e can  be  calculated  from  the  spacing  of  the  grids  and  the  dimensions 
of  the  grid  wires.  In  practice,  e must  be  adjusted  to  fit  the  characteristics  of  the 
instrument.  For  this  instrument  with  its  double  grids,  theory  gives  e < 0.  06,  but 
in  practice  e may  be  as  great  as  0.  12. 


For  the  DMSP  sensor,  let  us  assume  that  the  thickness  of  the  sheath  outside 
the  aperture  is  approximately  equal  to  the  grid  spacing.  Thus  the  applied  voltages 
of  interest  are 

' = 0 is  the  plasma  V = 0 

o 

i = 1 is  the  aperture  grid  Vj  = <t> , 

i = 2 is  the  retarding  grid  Vg  = V + 4> . 

i = 3 is  the  suppressor  grid  Vg  = V + 0 - 30 

where  V is  sweep  voltage  applied  to  retarding  grid. 

From  Eq.  (48),  the  effective  potential  on  the  aperture  grid  is 

<bA  = <t>s  (l  - i)  + t V (49) 

and  the  effective  potential  on  the  retarding  grid  is 

= <t>3  + V (1  - e)  - 30  e (50) 

The  flow  of  ions  to  the  collector  plate  is  regulated  by  the  grid  with  the  most 
positive  electrostatic  potential.  Table  2 shows  the  effective  potentials  as  a function 
of  the  applied  retarding  potential  for  <t>^  = -2  V and  e = 0.  06  and  e = 0.  12.  When 
<t>A  £ the  effective  potential  on  the  aperture  controls  the  current,  and  when 

^ the  effective  potential  on  the  retarding  grid  controls  the  current.  In  the 

data  analysis,  the  effective  potential  controlling  the  ion  current  should  be  used 
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Table  2.  Effective  Electrostatic  Potential  for  Planar  Sensor 


instead  of  the  applied  voltage  on  the  retarding  grid.  Use  of  the  uncorrected  retard- 
ing  grid  voltage  yields  ion  temperatures  that  are  too  high  by  10  to  40  percent. 


.i 


( 


,1 

. 


4.  PLASMA  SCALE  HEIGHT  DETERMINATION 

There  is  sufficient  information  from  the  analysis  of  the  electron  Mode  2 and  ion 
3 

Mode  2 data  to  compute  the  plasma  scale  height  Hp  at  a time  centered  on  the  two 

swept  voltage  modes  which  occur  17  sec  apart;  that  is,  Hp  is  determined  at 

(t  + 8.  5),  where  t is  center  of  an  electron  sweep,  to  be 
o ' o 


H (t  + 8.  5)  = k(T  +T)/M.m  g km 

P O 6 l P 


o 

where  g = acceleration  due  to  gravity  at  835  km  altitude  (km /sec  ); 

g = 7.  675  10'3  km/sec2 
Substitute  in  constants 


H (tQ  + 8.  5)  = 1.  083 


T + T 
e 

M. 


km 


(51) 


(52) 


3.  TELEMETRY  ALLOCATIONS  AND  ELECTRONICS 

With  the  exception  of  two  housekeeping  data  words,  all  the  data  from  the  SSIE 
experiment  is  transmitted  via  the  spacecraft  digital  data  systems  (OLS)  using  a 
single  NRZ  (nonreturn  to  zero)  data  signal.  This  signal  is  transferred  once  per 
second  in  phase  with  the  OLS  supplied  bit  clock  in  bursts  of  180  contiguous  bits  at 
a rate  of  1000  ± 1 bps  as  shown  in  Figure  11.  Upon  receipt  of  a read  pulse 
(SSIERED)  from  the  OLS,  the  SSIE  experiment  provides  180  data  bits  with  the  least 
significant  bit  occurring  first  in  the  first  word.  These  180  data  bits  form  twenty 
9-bit  data  words. 

The  outputs  from  the  three  experiments  are  equally  spaced  throughout  the 
1-sec  interval  and  arranged  in  the  twenty  data  words  as  follows; 


Type 

No.  Samples 

Words 

Electron  sensor  data 

7 

1, 

4. 

7,  10,  13, 

16, 

19 

Ion  sensor  data 

7 

2, 

5. 

8.  11,  14, 

17, 

20 

Event  monitor 

6 

3, 

6, 

9,  12,  15, 

18 

3. 


Rishbeth,  H. , and  Garriott,  O.  K. 
Academic  Press,  New  York,  p. 


(1969)  Introduction  to  Ionospheric  Phvsics, 
143.  ‘ 
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SSIERED  SSIERED  SSIERED 

Figure  11.  SSIE  Telemetry  Word  Transfer 

The  electron  and  ion  sensor  data  words  are  converted  from  a telemetry  output 

I voltage  to  an  amplifier  current  using  preflight  and  in-flight  calibrations.  This 

procedure  is  described  in  Section  5.  2. 

The  event  monitor  is  an  analog  indicator  of  the  voltage  applied  to  both  the 
electron  and  ion  sensors  during  their  respective  sweep  periods.  The  actual  volt- 
age to  the  sensor  components  is  derived  by  fitting  a straight  line  to  the  event 
monitor  analog  voltage  and  then  comparing  it  with  a preflight  calibration  of  applied 
volts  vs  analog  volts. 

The  event  monitor  also  indicates  calibration  sequences  and  "flags"  the  sweep 
periods  2 sec  prior  to  the  sweep,  as  an  aid  to  data  processing.  It  should  be  noted 
that  all  events,  such  as  calibrates,  sweeps  and  flag  voltages  first  appear  at  word 
18  of  the  20-word,  180-bit  burst  of  data.  Therefore  the  first  words  of  data  from  a 
sweep  or  calibrate  are  words  19  and  20  for  the  electron  and  ion  sensors, 
respectively. 

The  two  housekeeping  quantities  which  are  transmitted  using  the  spacecraft 
analog  system  are: 

i Word  48  - Output  voltage  of  thermistor  in  main  electronics  package. 

Calibration  given  in  Figure  12. 

Word  49  - Monitor  of  bias  voltage  applied  to  grid  of  electron  sensor. 

' Values  given  in  Table  1. 

<| 

5.1  Time  Sequence  of  Event* 

Figure  13  shows  the  sequence  of  all  functions  of  the  Event  Monitor,  Electron 
Sensor  and  Ion  Sensor.  At  all  other  times,  the  sensors  themselves  are  operating 
and  recording  data  in  the  fixed  voltage  mode.  Additional  information  is  listed  in 
j Table  1,  including  in-flight  calibrate  details  and  the  designation  of  the  telemetry 

output  for  the  electron  grid  resting  bias  voltage. 

1 

] 
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TEMP-*C 

Figure  12.  Thermistor  Calibration 
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ION 

SWEPT 

GRID 


I — / — — I — J 


time  <n  sec 


EVENT 

MONITOR 


■H  JLU 

4 V , 1 ' 


t VENT  MONITOR  PERIOD  A" 


NOTES 

1 ELEC  SENSOR,  DURING  BIAS  MODE  I,  + 1 5 V 
APPLIED  TO  OUTER  GRID  DURING  BIAS  MOOE  2, 
OV  APPLIED  TO  OUTER  GRID 

2 ELEC  SENSOR  COLLECTOR  IS  BIASED  AT  + 20  V 
WRT  OUTER  GRID  AT  ALL  TIMES  INCLUDING 
SWEEP  PERIODS 

3 ION  SENSOR  SUPPRESSOR  GRID  AND  COLLECTOR 
ARE  BIASED  AT  - 30  V AND  -20  V RESPECTIVELY 
WRT  SWEPT  GRID  AT  ALL  TIMES  INCLUDING 
SWEEP  PERIODS 

4 ION  SENSOR  APERTURE  GRID  IS  AT  SPACECRAFT 
GROUND  AT  ALL  TIMES 


SSI/E  TIMING  SCHEME 
DMSP  BLOCK  50  F2  SPACECRAFT  ONLY 

Figure  13.  SSIE  Time  Sequence  of  Events 


5.2  Preflight  Calibration  and  In-flight  Update 


Before  integration  with  the  spacecraft,  the  amplifiers  were  accurately  cali- 
brated for  response  to  known  input  currents  over  a wide  range  of  temperatures. 
The  preflight  calibrations  are  shown  in  Figures  14  and  15  and  in  Tables  3 and  4, 

for  a temperature  of  +20°C.  Both  amplifiers  have  a logarithmic  response  with 

-10  -5 

outputs  0 to  5 V for  ranges  10  -10  A in  the  electron  (negative  current)  ampli- 

-12  -7 

fier  case,  and  5 X 10  -5  X 10  A in  the  case  of  the  ion  (positive  current) 

amplifier. 

These  amplifier  responses  may  be  expressed  as 

Log(I)  = A + BV  (53) 

where  Log(I)  is  Log10  of  the  input  current;  V is  output  voltage  of  amplifier  to  tele- 
metry; and  A and  B are  constants. 

Using  e and  i as  subscripts  to  denote  the  electron  and  ion  current  amplifiers, 

from  the  preflight  calibrations,  we  obtain  the  following: 

A =-10.00  B =+1.000 
e e 

At  = -11.301  Bt  = +1.000 

Since  the  in-flight  data  received  can  commence  at  any  point  in  the  1024 -sec 
time  sequence  between  calibrates,  these  preflight  values  of  A and  B are  used  to 
process  the  initial  data  before  the  first  in-flight  calibrate  is  encountered  in  each 
data  set. 

As  shown  in  Table  1,  each  amplifier  is  fed  two  levels  of  known,  stable  input 
currents  once  every  1024  sec.  On  each  amplifier  these  consist  of  two,  2-sec 
periods  (14  data  points  each)  for  each  of  which  an  average  value  of  output  voltage 
V is  obtained. 

Let  Vle,  V2e,  Vj.,  and  V2i  be  the  calibrate  output  voltage  levels  for  Calibra- 
tions 1 and  2 of  the  electron  and  ion  sensors,  respectively.  We  now  obtain  updated 


values  for  A and  B in  Eq.  (53) 

Ae  - (6  V2e  - 9 Vu)/(Vle  - V2e)  (54) 

Be  = 3/<Vle  - V2e>  <55> 

At  = (7.301  V2l  - 10.301  - V2i>  (56) 

Bi  = 3/(Vli  " Vai*  <57) 


These  updated  values  of  A and  B are  then  used  for  all  sensor  data  in  both  fixed 
voltage  and  sweep  voltage  modes  of  operation  until  the  next  calibrate  is  evaluated. 
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Figure  14.  Electron  Amplifier  Laboratory  Calibration,  +20°C 


Figure  15.  Ion  Amplifier  Laboratory  Calibration,  +20°C 


; 


Tabled.  Electron  Sensor  Amplifier  I.  aboratorv  Calibration,  +20°C,  Block  50,  F2 


Calibration  1 (10  *'  A):  +3.964  V 
Calibration  2 (10‘9  A):  +0.978  V 

Monitor  Levels:  Electron  Sweep  +5  V -*  +1  V 

Sweep  Flag  +0.  6 V 

Calibration  1 +0.2  V 

Calibration  2 +0.4  V 


Table  4.  Ion  Sensor  Amplifier  Laboratory  Calibration,  +20°C,  Block  5D,  F2 

Calibration  1 (5  X 10'8  A):  +3.  984  V 
Calibration  2 (5  X 10  A):  *-0.076  V 

Monitor  Levels:  Ion  Sweep  +1  V — +5  V 
Sweep  Flag  +0.  8 V 
Calibration  1 +0.  2 V 

Calibration  2 +0.4  V 


